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SECÇÃO DE PROPOSTAS 


Em fodos os estabelecimentos SKF 
existem fécnicos e engenheiros espe- 
clalizados ao dispor de quem os queira 
consultar sobre tudo quanto diga res- 
peito a rolamentos, bem como a pres- 
tarem toda a colaboração na escolha 
dos mesmos e controlar a sua monta- 
gem. 
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Também, para um maior desenvolvi- 
mento do seu serviço fécnico, faz a 
SKF' um constante intercâmbio entre 
os seus estabelecimentos de todo qo 
mundo. 


SOCIEDADE SKF' LIMITADA 
LISBOA, Praça da Alegria, 66-A 
PORTO, Avenida dos Aliados, 152 
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TEODOLITOS — TAQUEÓMETROS AUTO-REDUCTORES — NÍVEIS 
DE ÓCULO—ALIDADES —BINÓCULOS — MATERIAL DE DESENHO 


INSTRUMENTOS DE FAMA MUNDIAL 
PELA PERFEIÇÃO TÉCNICA E PRECISÃO 


KERN & CIE, S. A., AARAU, SUÍÇA 


Representante em Portugal: 


EMÍLIO DE AZEVEDO CAMPOS & CA, LDA. 


CASA FUNDADA EM 1854 
P. Santo António, 137 - 1445 — Perto — Tel. 20254/5/6 
PR. Antero Quental, 47-4.º — Lisboa — Te). 553366 


PEÇA-NOS CATÁLOGOS DESCRITIVOS 


TÉCNICA — KXHIT 


C. D. U, 622.755 


Verificação da lei de Tillé para a distribuição 


dos caudais nos hidrociclones 


Pelo Eng.º de Minas (1.5. T.) NUNO GONÇALO BACHA DE ALMEIDA RIBEIRO 


1 — Introdução 


1.1 — Nos últimcs tempos tem-se procurado 
utilizar forças diferentes das da gravidade para 
provocar a hidroconcentração. É o que acontece 
nos hidrociclones. 

Sabe-se que, num fluido com partículas sólidas 
em suspensão, as forças gravitacionais provocam 
descida mais rápida das partículas de maior cali- 
bre ou das que, tendo o mesmo calibre e forma, 
tenham maior densidade, Uma maneira de activar 
este processo é fazer actuar sobre as partículas 
sólidas outras forças mais intensas que as da 
gravidade. Consegue-se este objectivo injectando, 
por meio de uma bomba, o líquido, com o sólido 
em suspensão, através de um orifício de entrada, 
tangencialmente à superfície interior do corpo 
cilíndrico ou cilíndrico-cónico do aparelho de 
concentração. É este o princípio do hidrociclone, 
no qual as partículas sólidas são, assim, subme- 
tidas a uma força centrífuga, cuja acção provoca 
hidroclassificação muito acentuada. 

Embora sejam completamente conhecidos os 
fenómenos que ocorrem quando uma partícula 
material isolada evolui num fluido homogéneo, 
infinito e isótropo, o mesmo não acontece se se 
trata dum meio heterogéneo, de numerosas par- 
tículas, confinadas num volume limitado, como é 
o caso do hidrociclone, pois aí teremos de contar, 
não só com a acção mútua entre as partículas, 
como com os efeitos de parede. 

Neste caso, para a determinação do compor- 
tamento das partículas há necessidade de recorrer 
a ensaios práticos, dada a grande dificuldade 
teórica de definir todos os parâmetros da ciclona- 
gem e de equacionar as leis quantitativas que 
regem os fenómenos postos em jogo. 


1.2 — Dado o interesse deste problema mon- 
támos, na Lavaria Piloto, do Laboratório de Pre- 


Assistente do 1.5. T. 
Ex-Bolseiro do 1. A.C. na €C, E. E.N. — Laboratório 
de Tratamento de Minérios 


paração de Minérios da C.E.E.N., um circuito 
contínuo para o estudo da ciclonagem, à escala 
industrial, de minérios portugueses de urânio. 

Assim, para perfeita conexão dos diferentes 
órgãos de tal circuito, julgou-se necessário deter- 
minar, em primeiro lugar, a influência de vários 
parâmetros no funcionamento do ciclone exis- 
tente, para que se conhecesse a respectiva distri- 
buição de caudais. 


1.3 — Embora exista uma considerável biblio- 
grafia sobre a utilização de ciclones para espes- 
samento e clarificação de polpas, o seu uso como 
processo de concentração de minérios é recente, 
nomeadamente quanto a minérios de urânio, 
onde a sua aplicação teve início em 1957, nos 
E. U. A. Infelizmente, não foram publicados rela- 
tórios com os resultados daquela aplicação. 

Os primeiros trabalhos importantes sobre ci- 
clonagem foram realizados a partir da última 
guerra, pela Dutch State Mines, em Limbury, 
Holanda (1), Esta empresa, nas suas minas de car- 
vão, utiliza ciclones de 350 mm para recuperação 
de areias de calibres inferiores a 50 mm, areias 
essas que são usadas como adensadoras do 
«meio» nas lavagens de carvão. 

No entanto, foi R. Tillé que deu a primeira 
contribuição para o estudo teórico do ciclone com 
a sua lei relacionando os caudais com os diá- 
metros dos orifícios de saída. 

M. Stass, no Congresso Internacional da Con- 
centração de Minérios, realizado em Stockholm 
em 1957, apresentou um interessante trabalho (1) 
procurando verificar as leis propostas por R. 
Tillé, onde conclui que a referida lei deveria ser 


(1) M. Stass (1957). The influence of the orifices on 
the Washing Characteristics of the Hydrocyclone. Inter- 
national Mineral dressing congress, Stockholm, 1958. 
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afectada de constantes que provavelmente de- 
pendiam das características geométricas de cada 
ciclone. Pareceu-nos assim lógico verificar, para 
o ciclone de que dispúnhamos, as leis propostas 
por R. Tillé e M. Stass e também determinar, 
para o caso concreto em estudo as leis que rela- 
cionam os diâmetros com as densidades do 
«meio» utilizado, do espessado e do transbordo 
do ciclone. 


2 — Montagem e técnica operatória 


2.1 — A montagem do ciclone-se fez conforme 
o esquema da fig. 1. 


Fig. 1 — Esquema da instalação experimental 
do microciclone de 1! 


A um depósito B ligou-se uma bomba tipo 
MONO, de 1 CV, destinada a impulsionar o pro- 
duto a concentrar para o ciclone A. Comple- 
tou-se o sistema com um manómetro F e com 
as válvulas de comando E e D, esta última inter- 
calada num by-pass. 


2.2 — Para a preparação de cada um dos en- 
saios procedeu-se do seguinte modo: 


a) Fechou-se a válvula E e abriu-se a vál- 
vula D. 


b) Introduziu-se em B água e o minério. 
c) Arrancou-se com a bomba €. 
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Bomba MONO-1CYV 


Ficava-se então em condições de iniciar os en- 
saios própriamente ditos (1). 


2.3 — Para a realização daqueles, abriu-se a vál- 
vula E e foi-se fechando a válvula D até se obter 
em F a pressão de trabalho desejada. Colhe- 
ram-se, em cada ensaio, duas amostras, em a e b. 


3 — Ciclone 


Na fig. 2 apresenta-se o desenho do ciclone 
utilizado nas experiências. 

Consta essencialmente de dois corpos, um ci- 
líndrico que se denomina «cabeça», e outro 
cilindro-cónico, que se denomina «corpo do 
ciclone». 

A cabeça do ciclone dispõe, superiormente, 
de uma tampa roscada e, inferiormente, é tapada 
por um septo metálico que tem ao centro um 
orifício (2) que se prolonga por um tubo. É pos- 
sível fazer variar o diâmetro deste orifício (2), 
introduzindo no tubo pequenas peças metálicas 
com orifícios calibrados. 

Na cabeça do ciclone há ainda um orifício 
para a saída dos produtos ciclonados, que se 
indica com o número (3). 

No corpo do ciclone, na sua zona cilíndrica, 
há um orifício de entrada dos produtos a ciclo- 
nar (4). O eixo deste orifício é tangencial à pa- 
rede interior do corpo do ciclone. 

No vértice do corpo do ciclone, dispõe-se de 
uma peça roscada onde podem ser adaptados 
bicos com os diâmetros de saída convenientes. 

O orifício de saída do ciclone é indicado pelo 
número (1). 

Daqui em diante, tudo o que se refira aos dife- 
rentes orifícios do ciclone será afectado com os 
convenientes dígitos 1, 2, 3 ou 4. 


4 — Numenclatura usada 


E conveniente definir os diferentes símbolos 
usados. 
Assim temos: 


1 — indica o orifício de descarga dos es- 
pessados 
2 — indica o orifício de descarga do trans- 


bordo 


(1) — Se se pretendesse ensaiar com «meio denso», 
criado por exemplo por magnetite, para aumentar a den- 
sidade da polpa, a introdução desta seria feita em B, com 
a bomba C em funcionamento, para se obter um meio tão 
homogéneo quanto possível, 


o] 
— 
mm 


ANY 
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Fig. 2 — Microciclone de 1! 


Po — pressão de trabalho do ciclone 

D; — diâmetro do orifício de descarga dos 
espessados 

D: — diâmetro do orifício de descarga do 
transbordo 

Qo — caudal total da alimentação 

OQ, — caudal do espessado 

OQ» — caudal do transbordo 

do — densidade da alimentação 

di — densidade no orifício de descarga 

ds — densidade do transbordo. 


S — Condições de realização dos ensaios 


Tillé e Stass verificaram que a densidade do 
fluido utilizado na ciclonagem tinha pequena 
influência na lei que relaciona D;/D: com Q,/0s. 
Assim, nas experiências realizadas, utilizou-se 
unicamente água. 


Os diâmetros dos orifícios fizeram-se variar, 
utilizando bicos com os diâmetros seguintes: 


A pressão de alimentação foi de Po == 2,5 
kg.cm”?. 


6 — Resultados 


Os resultados obtidos constam do quadro 
na página seguinte: 


7 — Análise dos resultados 


7.1 — Com os valores obtidos traçou-se, em 
coordenadas logarítmicas, a curva que consta da 
fig. 3. 
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Fig 3 — Variação de D;/D; em função do Q,/Q» 


Esta curva, entre os valores de Di/D: 0,65 e 
1,4 pode considerar-se uma recta, pelo que é 
possível escrever 


ou 
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2,000 
2,480 


0,0304 
0,0380 
0,0722 


0,113 
0,175 
0,182 
0,274 
0,311 
0,497 
0,645 
1,000 
1,338 
1,610 
2,294 
3,880 
4,070 
7,400 
8,680 


17,500 
19,347 


| 2,960 21,929 | 


R. Tillé propunha, para aquela zona da curva, 


a fórmula 
OQ, D an | 
Qs ( a. Ds (1) 


indicando para a e n os valores, respectivamente, 
0,921 e 3,5. 

Para o nosso ciclone, conforme as determina-. 
ções que efectuámos, obtivemos 


n== 3,757 
a == 0,990 


que, substituídos em (1), dão: 


O ( D: (2) 
Q: 0,990. Ds, 

As deflexões que se observam para cada lado. 
do intervalo compreendido entre os valores de 
D:/D: 0,65 e 1,4 corresponderam a ter-se, para 
pequenos e grandes valores de Di / Ds», verificado, 
num dos orifícios, a entrada de ar em vez de 
saída de polpa. 


Notou-se a ocorrência de dois casos: 


a) Com Dy muito pequeno e Ds muito largo. 
todo o caudal saía por Ds, verificando-se 
sucção de ar através de Di. 


b) Com D: muito pequeno e Di; muito largo. 
verificou-se o contrário, com toda a saída. 
por Di. 


Entre estes dois casos, o caudal é distribuído 
entre os dois orifícios de acordo com a lei ex- 
pressa em (2). 


7.2 — Sendo: 
Q: + Q: ss Q, 


O: di + Q: d: = Q, do 


ou 


di == do + (do — ds) E 


pode escrever-se, de acordo com (1), 


dida fo = a ( o ) 


ou seja 


Di 
Ds» 


1 
do—d: ) 3,757 (3) 


== 0,990 ( 
di — do 


8 — Conclusões 


Pela fórmula (3) atrás deduzida é possível, 
para cada valor d, da densidade da polpa a ciclo- 
nar, determinar a curva de variação da relação 
D:;/D: dos diâmetros dos orifícios de descarga 
do espessado e do transbordo do ciclone estu- 
dado com a densidade dy do produto espessado. 

Torna-se assim possível seleccionar os pares 
valores (Di /Ds) que permitem obter, com aquele 
ciclone, uma densidade de espessado dy corres- 
pondente à melhor recuperação. 
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NOTAS INFORMATIVAS | C. D. U. 624.311,5/94 


Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 


— Elementos extraídos das estatísticas mensais do Repartidor Nacional de Cargas ( R.N.C.) 


Nora: As produções e os consumos das empresas do R. N. C. representam 
cerca de 94 º/, dos totais do Pais. 


JANEIRO 


I — Breve nota mensal 


Sob o ponto de vista hidrológico. o mês de Janeiro 
apresentou-se, no conjunto, muito acima da média. 


[HI — Elementos gerais (GWh) 
a) Mensais 


| Variação 
0/y 


[Em 162 


| 
Produção hidráulica (Ph)... 502,6 | 322,3 | 4 7 


Produção térmica (Pr) ..... 0,0 0,0. U 
Produção total (PT)... ...| 902,60/ 3223 + | 
Exportações (Ex) . cc... .. = — 
Importações (1) « . cw «sue o — | 0,0. — 
Produção para con- | 
sumos não perman. (Penp). 554 


C0,3 | + 9 
Produção para con- | 
sumos permanentes (Pcp) (!) | pa 292,0 | + 6,0(' | 
ad a IV — Energia armazenada nas principais albufeiras 


NOTAS: 
(1) A produção para consumos permanentes (Pcp! é determinada | | No fim do mês 
pela seguinte expressão: Pcp = Pr — Penp + | — Ex Albufeiras ; A 


| GWh | “o (1) 
(2) O aumento percentual de produção para consumos permanen- a | 


tes, tendo em conta a incidência dos domingos e dias especiais, é | 


de 78 9o- Paradela. . . . cc. ..«| 2204 99,1 
Venda Nova . .. «+» «| 1026 80,2 
III — Diagramas de carga dos dias característicos Salamonde . . +... EUA a 23,5 85,1 
Canigada cs usviccic] 819 95,8 
4.º feira; CONES ses usas 5 q | 8143 94,6 
18-1-961 | 17-1-962 Sestolo do M6d. wo +» | 151,8 93,1 
ai air es ns Gullhofrel «cw canso 7,5 90,4 
Produção hidráulica (Ph) — MWh 10480 | 11367 Lagoa Comprida * Cn & à | 29,90) | 81,6 
Produção térmica ( Pr)— NWh.. () | 0 Bata Euald = ss vis sá | 58,2 91,2 
Produção total (PT) — MWh . ..| 10480 | 11367 Pracana ME O DL E A A E DO 11,4 88,4 
Trocas com [ Exportações MWh . — — PÓvOa «verve a 8,9 (*) 2,1 
Espanha | Importações MWh . — — Total. . .| 953,8 91,4 
Produção para consumos não per- | 
manentes — MWh. . . ..... 1644 1454 
Potência máxima MW ....... 531,5 | 638,5 Notas : 
Potência mínima MW ....... 241,0) | 2680 | (1) Coeficiente de enchimento em energia das albufeiras, 
ra É ; do 
cap nha ço a | dá (2) Ineloi 1, GWh armazenados em Vale do Rossim no 
À | fim do mês, 
ido e, Ela coesa! 0,41 0,42 , VW Ei 
Pot. máx. | E ' (*) Inclui 1,4 GWh armazenados no açude do Poio no 


fim do mês, 
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VALVULAS 


PASSAGEM e RETENÇÃO e REGULAÇÃO e CUNHA 


PARA LÍQUIDOS, VAPORES 
E GASES INCL. NHs E SO» 


DN 10m/m ATÉ 400 M/m 
PN 6 a 40 KGS / cuz 


CONFORME NORMAS DIN E ASA 


EXECUÇÃO EM FERRO FUNDIDO, 

FERRO FUNDIDO C/GRAFITE 

ESFEROIDAL, AÇO VASADO 
E AÇO INOXIDÁVEL 


VÁLVULAS DE SEGURANÇA, DE 


REDUÇÃO, COM DIAFRAGMA. ELECTRO- 
MAGNETICAS E VÁLVULAS ESPECIAIS. 


ESTUDOS: 
'PRUJEGTOS: 
ORÇAMENTOS: 


SOCIEDADE PORTUGUESA DE VÁLVULAS, 


RUA ACADEMIA DE CIÊNCIAS, 5 - LISBOA 
TCETEPONES: IAC PE IIZAB- TEST CL 


PORTO: PRAÇA D. JOÃO |, 25-10 - TELEFONE 24771 


VALVULAS DE QUALIDADE PARA TODAS AS APLICAÇÕES 
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TRANSFORMADORES 


LICENÇA SIEMENS SCHUCKERT 
ATÉ 1600 kVA 35 kY 


MOTORES TRIFÁSICOS ASSINCRONOS 
ATÉ 120 Cv 


EMPRESA NACIONAL DE APARELHAGEM ELÉCTRICA, 3. À R. Lav. 24 DE JULHO, 158 — LISBOA 
oistrisuiDORES 7 CSA ES dA EQUIPAMENTOS ELÉCTRICOS, S. A. R. L. 
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C.D.U. 624.39;517.56 


DETERMINAÇÃO DO ESPECTRO DE POTÊNCIA 
A PARTIR DE REGISTOS CRONOLÓGICOS 


FERNANDO EDUARDO REBELO SIMÕES 


Centro de Estudos de Electrónica da Comissão 
de Estudos de Energia Nuclear do Instituto de Alta Cultura 


Assistente do Instituto Superior Técnico 
SUMÁRIO 


Faz-se o estudo comparativo do método de análise directa de frequência e do método de autocorrelação, para a determinação do espec- 
tro de potência a partir de registos cronológicos, mostrando-se cômo a própria limitação dos registos conduz naturalmente à noção de 
«anela» utilizada nos referidos métodos por Blackman & Tukey (1), e por Ackasu (2). 


SUMMARY 
The autocorrelation analysis and direct Fourier analysis are compared, for the determination of the power spectrum from time records. 
lt is shown how the limitations of the records give rise to the concept of «window» which is used in those methods by Blackman & Tukey (1), 
and by Ackasu (2). 


Por se tratar de um assunto actual e de grande interesse tem-se ocupado o Centro de 
Electrónica, desde há algum tempo, com os problemas da Análise Harmónica Generalizada. 

A aplicação e o interesse que estas técnicas têm encontrado, nomeadamente por parte 
do Laboratório Nacional de Engenharia Civil, têm constituído um forte estímulo para o 
prosseguimento deste estudo. 

O presente artigo dá continuidade aos artigos publicados pelo Prof. Abreu Faro (º) 
sobre o mesmo assunto, e resulta justamente da necessidade de precisar de modo quantitativo 
os métodos de análise preconizados nas referidas publicações. 

Este trabalho integra-se num projecto de investigação científica sobre Análise Harmó- 
nica Generalizada e nele se utiliza material de laboratório que foi adquirido através de um 
subsidio que a NATO concedeu ao Laboratório de Telecomunicações do 1. S. T. para este 
projecto de investigação. Além deste material utiliza-se ainda material pertencente ao Centro 
de Estudos de Electrónica. 


1. INTRODUÇÃO 


1.1. Um conceito fundamental na teoria dos processos estocásticos, e de grande interesse nas 
aplicações, é o de densidade espectral de potência, ou como é mais correntemente designado, espectro de 


potência. 
Esta grandeza, W (f), pode definir-se duma maneira estatística, a partir das funções si (t) 
representativas do processo, e supostas de duração ilimitada. 
T | 


2" 2 
A expressão do conteúdo quadrático médio (no tempo) de s/ (t), relativo ao mesmo intervalo T) 
será dada por 


Designemos por si. (t) a função si (t) restringida ao intervalo (— 


1 T/2 
fas Í [si (DJ? dt (1.1) 
TJ rp 
(1) Ref. 11] 
(?) Ref. [2] 
(*) Ref. [3] 
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Consideremos agora a nova função s'i (t), definida em todo o domínio: 


tod) gcT 
2, 
st (t) = (1.2) 
| Ps 
O 2 


e seja Si. (f) o seu espectro de Fourier no domínio da frequência. 
Utilizando a relação de Parseval, podemos escrever 


Pi=[ Wi(g de (1 


 —DO 
em que 
1 3 
wi =| 54 | 
ou, fazendo intervir os dados, 
à LAR o, =jemte)|* 
Wa (1) = 7 | sr (t) e dt (1.3) 
— T'/2 
1.2. Consideremos agora as duas operações seguintes 
— passagem ao limite, quando T tende para infinito ...... lim 
T -—» oo 
— média, no conjunto das funções s' (t) ......cccc.... a 
Verifica-se que se tem em geral para os processos estocásticos as seguintes relações: 
< lim Pi> = lim jo e == P 
<limwr(0> &lim<WwO>=W(f 
T-—+ 00 T-+00 
por W. (f) não convergir geralmente para um limite. Pode finalmente escrever-se 
P = W (£) df (1.4) 
À 1 TR —ij2% 2 
W= lim — < ) s4 (t) e nz “dt > (1.5) 
T-o T —T| 


em que P representa a potência média, no tempo e no conjunto, do processo, e W (f), a densidade 
espectral correspondente, isto é, aquilo que se designa correntemente por espectro de potência e 
que nos propomos calcular. 


1.3. Um outro conceito igualmente importante é o de autocorrelação, que se pode definir 
para cada função s! (t) como 


T;2 / 


Wi (5) — lim es (=), t+) 1.6 
Teo T) a dd (1.6) 


(1) Em geral nos processos físicos pode fazer-se corresponder |s (t)]? a uma potência. Daí resulta que Pr àparte 
um factor terá o significado duma potência média e Wr (f) representará portanto uma densidade de potência média, 
no domínio da frequência, de — os a + 05. 
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O teorema de Wiener-Khintchin afirma que a média no conjunto das autocorrelações é a trans- 


formada do espectro de potência atrás definido. 
jamtz 


YO=<rig>= [" was” q (1.7) 


jamfr 


W (£) = ) Gm ST radar (1.8) 


Nos processos ergódicos, nos quais as médias tomadas ao longo do tempo coincidem com as 
médias no conjunto, as funções 1º (7) são todas iguais: a autocorrelação '” (7) pode ser calculada 
a partir duma única função, e o espectro W (f) pode dela ser obtido, através de (1.8), sem implicar 
uma média no conjunto. 

Esta será necessária, porém, na determinação directa (1.5). 


1.4. Na análise que vai seguir-se suporemos sempre os processos estacionários e ergódicos, 
definindo-se portanto a autocorrelação a partir de (1.6), mas suprimindo o índice i por desnecessário (1). 


Será, assim, ” 
T/ 
F (5) = im -— (1 sn) .s (+ + e dt (1.9) 
T+ oo - —T/2 2 


Esta função pode ser definida estatisticamente, em virtude da ergodicidade admitida, 


FP()=<s(t.sd(t+)> (1.10) 


sendo independente do instante +. 


2. AS LIMITAÇÕES DAS ESTIMATIVAS DE W(f) E W(z) 


2.1. Do que se disse acima resulta que a determinação do especto de potência a partir 
de registos cronológicos — nosso objectivo — se encontra limitada naturalmente pela impossibilidade 
de se dispor de registos de duração ilimitada e de registos em número infinito. 

Esta limitação vai acarretar, como veremos, uma dupla imprecisão das nossas estimativas: 


— imprecisão devida ao facto de ser impossível colher informações exclusivamente sobre 
uma dada frequência f 
— imprecisão devida ao carácter aleatório que resulta para a estimativa. 


2.2. Para ilustrarmos estes factos, consideremos um registo si(t) contínuo, definido entre 
E = É ma : es sai 
os instantes — E e EL (Fig. 2-1) e procedamos ao cálculo do espectro de potência do processo de 


que si (t) faz parte, admitindo que T pode, em princípio, ser tão grande quanto quisermos. 


As expressões (1.5) e (1.9) apenas poderão tomar respectivamente a forma 


Ti — j2nft jo 

Wi (£) = si(t)e dt (2.1) 

—T 2 

(T—|T|)/2 
1 A TR T 
vi= o fa(i- Dei (e+ Da 

E [m 3)” 2 (2.2) 

—(T—|7/)2 

(1) Este será usado mais adiante, como veremos, para caracterizar determinadas estimativas de W (x). 
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Ta 


Fig. 2-1 


Estas expressões já não se podem obter uma da outra por transformada de Fourier (1), sendo 
a segunda apenas definida para |7| < T. 

Daqui resulta que se quisermos obter uma estimativa de W (f) através duma estimativa da 
autocorrelação, há que estender o domínio desta função para além do intervalo (— T, T), a fim de 
que se lhe possa aplicar o integral de Fourier. 


Duas alternativas se nos oferecem: 
Ed º T T A . 
— ou se estende s.i (t) para além do intervalo (— FL de) em que é definida 


— ou se estende Wi (7) para além do seu domínio (— T, T). 


2.2.1. Extensão do domínio de si (t) 


Façamos a extensão definida em (1.2). 
A expressão (2.1) não se altera, mas a sua transformada em 7 escreve-se 


nO =) (ts (t+) dt, com + (-00,00) (2.3) 


dela se podendo obter W.i (f), segundo a expressão 


o 
+ 


— j2nfr 
Wi(£) = / Wr()e ds: (2.4) 
O cálculo de Wi. (f) através de (2.1) afigura-se menos trabalhoso do que (2.4), e será em 


princípio preferível, a menos que se pretenda obter especificamente uma estimativa da autocorrelação. 
Esta estimativa está relacionada com “Vi (x), (2.2), em virtude de (1.2), tendo-se 


Wi (5) = Dr (5) HÁ (5) (2.5) 
PE TIS T 

Dr (7) = T (2.6) 
0 [> T 


Fig. 2-2 


(1) Notar que (2.1) coincide com (1.3), pelo que Bm < WA ()>=Wt(f) 
+ OO 


Também lim WÃi()=<vi()D="T(5) 


TOC 


TEONICA 
404 


| 
| 
| 
| 
| 
| 
em que | 


2.2.2. Extensão do domínio de WA (7) 


Tomemos a função 


E (6) tI<T 


pra (5) = 
Ú s|>T 


(2.7) 


A expressão (2.2) não se altera dentro do intervalo (— T, T), mas passa a anular-se fora dele, 


vindo a sua transformada em f 
oo 


' — j2 fz 
WHO = | Wit de (2.8) 


— (HO 


que agora não coincide com (2.1). 

Estas estimativas (2.1) e (2.8) estão porém relacionadas através dum integral de convolução (1), 
devido à relação (2.5). 

Com efeito, passando esta igualdade para o domínio da frequência e pondo-se 


ia -—jJ27fT 
Or (8) = | Dr (7) e d (2.9) 
resulta finalmente 


W+ (8) = Qr (8) + W'+ (f) inasei 


= [Qr(t-E)W'+ (F) dê 


Esta fórmula relaciona as duas estimativas de W (f) que se obtêm por transformada de Fou- 
rier das duas estimativas de '” (7) provenientes das extensões de domínio referidas. 


3. AS “JANELAS”. O PROBLEMA DA ESTABILIDADE 


3.1. Qualquer dos processos atrás descritos conduz a uma estimativa do espectro de potência 
que está sujeita a flutuações estatísticas que interessa reduzir na medida do possível. 

Evidentemente que esta redução se pode obter efectuando várias estimativas e tomando a sua 
média aritmética. É de resto a única solução no primeiro dos processos apontados. 

O segundo processo porém sugere-nos uma maneira diferente de aumentar a estabilidade 
da estimativa, por eliminação dos valores calculados de Wi (7) para os quais a instabilidade 
é maior. 

Trata-se, como é evidente, dos valores correspondentes aos atrasos | 7 | na vizinhança de T. 
Quer dizer que em vez de anularmos '[4. (+) apenas para |7| >> T, o devemos fazer para |7| > Tm, 
com Tm < T. 

Esta eliminação de elementos, sendo realizada num plano intermédio não é obviamente aplicável 
no primeiro processo, em que uma eliminação de dados só se poderá fazer em si (t), com 
a correspondente redução de informação à custa de T. E o valor finito de T constitui justamente 
o inconveniente que os métodos que estamos descrevendo procuram eliminar. 

Claro que comparando neste ponto do nosso raciocínio os dois métodos, é de esperar uma certa 


(1) Convolução entre Fix) e G (x): 


F()* Gl)= [Fix —y) G(y) dy= [F(y) G(x— y) dy 


— O — (O) 
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equivalência, uma vez que ao cálculo mais simples pelo processo directo (o 1º) se contrapõe a neces- 
sidade da realização de várias estimativas, ao passo que o processo indirecto (o de correlação), 
envolvendo embora uma transformação de Fourier, pode utilizar apenas uma estimativa, em igualdade 
de condições de estabilidade, como veremos adiante. 


3.2. O processo directo pode considerar-se como consistindo em se multiplicar a função 
si (t) — de que se conhece apenas si. (t) — por uma função B, (t) assim definida : 


B, (9 =| [A 1 (3.1) 


0 |t|> T/2 
função esta que se comporta como uma «janela» no domínio do tempo (1). (Fig. 3-1) 


Anâlogamente o processo da autocorrelação pode considerar-se como consistindo no produto 
de W14. (7) por uma «janela» no domínio dos atrasos (2) D, (+) definida como segue 


|t|<Tm=T 


(3.2) 
|t]|>Tm=T 
ofo) 
3 0 a 
Fig. 3-1 Fig. 3-2 


3.3. Podemos agora generalizar estes métodos, considerando em vez das janelas unitárias 
Bo (t) e D, (7), respectivamente as janelas B (t) e D (7), funções pares anulando-se para |t | > T/2 
e |7| >> Tm,(3) tal como as que acabamos de descrever, mas de resto inteiramente arbitrárias. 
(Figs. 3-3, 3-4). 


(o) 


m 


Fig. 3-3 Fig. 3-4 


4 


O critério que há-de presidir à sua escolha baseia-se, como veremos, em considerações de 
simplicidade, estabilidade, precisão e, no caso dos registos discretos, de largura de banda do processo. 


(1) Segundo a denominação de Blackman & Tukey: «data window». 
(2) idem: «lag window». 
(3) Tm &, T. 
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Vamos agora considerar separadamente cada um destes métodos, começando pelo método 
indirecto. 


4. O MÉTODO DA AUTOCORRELAÇÃO 


4.1. CASO CONTÍNUO 


4.1.1. Seja um registo contínuo de duração T, da função s! (t) (Fig. 2-1) pertencente a um 
processo estacionário e ergódico, e efectuemos a estimativa (2.2). 


(T—| = |)/2 
1 T T 
Wi ms = 2.— i PRESS) i Re 
no= a / sp (1 3 st(t+ 5) 
(E | 2] ja 


Multiplicando em seguida Wi (r) pela janela D'(7), de largura 2 Tm, como foi indicado, 


obtemos o valor corrigido 
Vi()=D() Vil) (cf2.4) (4.1) 


ao qual aplicamos a transformada de Fourier. 


OCO Ci 
— |2nf 7 
Wi(g)= f Vit) e dr (4.2) 


Este valor W'(f) é, tal como Wi e Wi, uma variável aleatória e constitui a nossa estimativa 
do espectro de potência W (f). 

Há pois que estudar W'(f) do ponto de vista estatístico, não só para relacionarmos esta 
grandeza com o espectro pretendido W (f), mas também para determinarmos as flutuações que 
a estimativa pode sofrer. 


4.1.2. Suponhamos então que era possível efectuar uma estimativa W'(f) por cada uma das 
funções s! (t). 
Teríamos, nestas condições, 


a — j2rfr 
<wú>= [o D()<TilD)> e dr, 
ou seja 
— j2xf: (2) 
<Wi(f)> -f D(7)WY(7) e dr, comr(—co,os) (4.3) 


em virtude de o processo ser suposto ergódico [v. (1.10) e nota ao $ 2.2.]. 
A expressão (4.3) representa a transformada do produto D (x) '” (+), ambos os factores admi- 
tindo transformada (o que não acontecia em (4.1). 


oo 
(1) Wi(f)=— af wi (7) cos 2x fr dr, por ser Wi(z) uma função par. 


o 
(2) Embora <p d. (=) > seja definido só para |+| < T, o factor D (7) permite estender 7 a todo o domínio. 
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Teremos então 


| <Wi()>=0Q(f) * W(£) (4.4) 
ou seja DES SAE 
< Wi (ti) > = fo (Ei —£) W (£) df (4.5) 
em que amo 
= — j2nfz 
Q()=|D(m)e dr (4.6) 


A expressão (4.5) mostra-nos que a média das estimativas de W (f) relativas a uma dada fre- 
quência fi não coincide exactamente com W (fi), mas resulta pelo contrário da contribuição de 
outras frequências, através da «janela espectral» O (f). 

Esta contribuição será tanto menor quanto mais «aguçada» for a forma de Q (f) (Fig. 4-1). 

Em particular, seria 


<Wi'i(t)>=W(f) para Q (= (f) 


Claro que esta hipótese nunca se verifica, dada a circunstância de D (7) ter uma duração 
limitada. 

Estes factos implicam a existência de imprecisão na estimativa (como já enunciáramos no $ 2.1) 
no sentido em que aquilo que se calcula, mesmo em valor médio, diz respeito não só à frequência 
pretendida, mas também a outras. 


4.1.3. Resolução 


Tem mais interesse para o projecto o conceito de resolução no domínio da frequência, conse- 
quência imediata da referida imprecisão, e que se pode medir pelo afastamento mínimo 4f de duas 
determinações Wº (fi), W'! (fi + 4Af), que não conduza a uma sobreposição apreciável das frequências 
em jogo nas duas estimativas. 

Para definir 4f, lembremo-nos que as estimativas W! (f) têm transformada no domínio de 7 


de duração limitada (|7| << Tm), bastando portanto conhecê-las nas frequências múltiplas de coli 
(teorema do sampling), pois ficarão assim completamente definidas. 2 Im 

Dado que porém as janelas Q (f) habituais introduzem para esta separação forte sobreposição 
de duas estimativas contíguas, será preferível tomarmos como relativamente independentes as esti- 
mativas correspondentes a uma separação dupla (fig. 4-2). 


Será pois este valor isa <A (4.7) 


/ 


que definiremos como resolução. 


we) 


À 


N 


Fig. 4-2 


| 
| 


)j 
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4.1.4. Alguns tipos de janelas D (7), Q (f) 


E ia 1 |7|<Tm ei 
O É de j 
(a) j Pet ia 
à sen 27f Tm co 3 Ta fa E iá 
O (f) E NLRÊ RA (4.9) 


Estas janelas, com Tm==T correspondem à atitude simplista descrita no 8 2.2.2 [cf. (3.2)]. 
Vemos (Fig. 4-3) como os lobos laterais apresentam amplitudes apreciáveis (=20º/ do lobo prin- 
cipal). 

ae T, 
Da (5) = | E (4.10) 


(b) 0 |T|> Tm 


(4.11) 


sen Tf Tm ) 


Fig. 4-3 


Estas janelas, com Tm = T correspondem à atitude descrita no $ 2.2.1 [cf. (2.6)] (2). 
OQ: (f) apresenta lobos laterais menos acentuados que Q, (1), mas tem o lobo principal mais largo. 


E Cd [7] < Tm 
Ds» (7) == 2 Tm 


(c) ' || > Tm 


(4.12) 


À 1 1 1 1 
oy=Lam+rto it) + 0 eo (4.13) 
E ã =“ 4º 2 Tu 4" pts 
(1) O valor de Do (7) para += + Tm não intervém em Qo (f) após a integração. Pode portanto ser qualquer 
finito. Notemos porém que quando Ds (7) se obtém de Q, (f), fica definido em 7==+ Tm, com o valor 1/2, 
(2) D, (7) equivale a Bo (t), (3.1), quando Tm = T. Trataremos desta equivalência no método directo. 
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Estas janelas, denominadas Hanning, apresentam lobos laterais mais baixos (= 2º/o do prin- 
cipal) e mais estreitos, podendo estas características acentuar-se ainda mais à custa duma maior 
complexidade das janelas usadas. 

Dum modo geral podemos dizer que Q» (f) representa um bom compromisso entre a necessi- 
dade de estreitar o lobo principal (o que exige D (7) com andamento semelhante a D,) e a neces- 
sidade de diminuir os lobos laterais, a fim de obviar a um possível acidentado do espectro (o que 
exige D (7) variando gradualmente). 


4.1.5. Estabilidade 


A estabilidade das estimativas pode medir-se pela sua dispersão relativa Ac, que se define como 
sendo o quociente entre o desvio quadrático médio de Wi (f) (variância) e o quadrado do valor médio. 
Estas médias são evidentemente tomadas no conjunto. 
Assim 
var [Wi (fi) ] (4.14) 
o = Re =— A 
<W' (fu) > 
Pode-se fazer uma ideia aproximada da distribuição das estimativas, comparando-a com uma 
distribuição em y? com k graus de liberdade, que apresenta a mesma variância e a mesma média. 
Como para estas distribuições se tem 


k As = 2, (4.15) 


conclui-se que podemos descrever a estabilidade através deste número k. 

Um conceito equivalente e certamente mais sugestivo é o de dispersão D, — em decibels — 
correspondente a uma dada fracção p da distribuição, e que pode ser obtida de tabelas das 
distribuições em 32. 

Por exemplo, as tabelas mostram que com k — 4 e p = 80 */ se tem Dp = 10 db. 

Quer dizer que neste caso a probabilidade de que uma estimativa esteja compreendida numa 
faixa de imprecisão, de 10 db entre valores extremos, é de 80 º/o (1). 

Existem expressões simples que dão a relação entre k e Dp de modo aproximado. 


Assim pode escrever-se “A 
k=1+ — (4.16) 
p 
tendo-se 
QUADRO I 


A | 35 100 | 250 | 400 | 625 | 840 


4.1.6. O cálculo de 47 pode fazer-se na hipótese de si (t) ter distribuição gaussiana. 
Blackman & Tukey (2) mostram que se podem obter, nesta hipótese, os seguintes resultados: 


ffiat-m+o err were we (ESPETO) arde 
o 2rr 1 
Ao (p= 0000000... — E, 
| fo (E— £) W (f) af] 
sendo T=T-eT.(a. o (4.18) 


(1) O valor médio fica mais próximo do extremo superior, tanto mais quanto menor for k. 
(2) Ref [1] — pág. 16 e seguintes. 
pág. 100 e seguintes, 
(3) « é um factor introduzido com o fim de tornar viável o cálculo de var [wi (1 , O qual desperdiça alguns 


dados na determinação de 4 A (7), e cujo valor indicado pela prática anda à volta de 1/3. 
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No caso de se usarem p registos de duração Tk, para cada um dos quais se calcula Wi (f), e 
de que se toma depois a média aritmética, dever-se-á pôr na fórmula anterior 


T'=p(Tk—a«Tn)=T-paTm (4,19) 


A expressão (4.17) é pouco prática, devido à sua complexidade, mas é susceptível de se sim- 
plificar consideravelmente, em dois casos: 


— Se o espectro W (f) variar lentamente, na vizinhança |1/T'| da frequência fi, tem-se 
aproximadamente, para as janelas Q (f) habituais, 


Tm 
Ao (fi) = 4.20) 
(f1) ( 
para fi > se , crescendo este valor à medida que fi diminui, até atingir o dobro para f=0 
m 
Falando em termos de graus de liberdade, será (v. 4.15) para este caso 
fia 2 ts E 
Tim Tm 
4.21 
- (4.21) 
k — fi > O 
m 


— Se W (f) for constituído exclusivamente por um pico de potência na frequência f, 


1 
(e em — f,, claro), de largura pequena comparada com en mostra-se que se tem neste caso 


Às == 1, ousejak = 2, para qualquer frequência fi. 


Este resultado verifica-se para W()=A[5(f-f)+o(f+Hf], o que constitui um 
caso extremo. 

Na realidade a existência dum pico sobreposto a um espectro acarreta uma dispersão relativa 
Az igual à unidade apenas na zona em que as janelas Q (fi + £,) tomam valores apreciáveis ; isto é, 


Tm 
Az ==1 na região do pico, mas tende para e à medida que nos afastamos deste e nos aproxi- 


mamos de zonas onde W (f) varie de modo mais gradual. 


4.1.7. Vemos assim que a estabilidade se opõe em parte à resolução. 
Para baixarmos esta há que aumentar o atraso máximo Tm. Este aumento, se não actuarmos 
na duração T, irá por sua vez aumentar a dispersão. 
curioso notar como estes conceitos obedecem a relações de compromisso, do tipo das relações 
de incerteza da Física. 
Com efeito, e atendendo a que para k suficientemente grande é T' = T, podemos sintetizar os 
resultados anteriores sob a forma 


Af, Ag = (4.22) 


a 
T 


1 
em que E é uma constante para cada registo. 


42 CASO DISCRETO 


4.2.1. Suponhamos agora que dispomos de n 4-1 dados, igualmente espaçados de dt, e que 
podem ser provenientes duma amostragem (sampling) efectuada sobre um registo contínuo (Fig. 4-4), 
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Tal como no caso contínuo, começamos por determinar um valor médio (no tempo) dos produ- 
tos cruzados 
H T T 
sT(t—— | sT ( + =) (4.23) 
2 2 


Para maior simplicidade vamos dispor os valores conhecidos de si (t) a partir de t=0, ten- 


do-se assim o conjunto de valores si (q At), com q=0,1,2,...n 


As nossas estimativas a (7) tomarão agora a forma 
1=0,1,2....n(1) 


0 —f 

à 3 si (qát) si (qát + rát) 
F n+1-—r q=o 
Er (rá t) = (4.24) 

1 E q 

—— >> s si (qAt) si (qiát+4rÃt) t=0,—1,— 2 coil 
14) 
T=nat 
t 
0 


Fig. 4-4 
P - P i * - -— 
sendo portanto "+ (r At) uma função par como o era VY (7), mas constituindo uma sucessão de 


valores discretos, igualmente espaçados de Át, 
4.2.2. Multipliquemos estas quantidades pelas janelas D (7), de largura 2m À t, sendom < n. 


Obtemos assim as estimativas corrigidas 
(4.25) 


Wi (rÃt) = D (5) Wi (r At) 


O passo seguinte seria achar a transformada de Fourier desta função. No caso presente porém 
estamos a tratar com funções que só não são nulas em intervalos de medida nula, pelo que vamos 


multiplicar os valores anteriores por funções à afectadas de coeficientes adequados. 


Obtemos assim 
a E ao; 
PU tag 5 Vi(mABI(G+EmA)+ 3 Wi(ADIG-rAD+ 
2 t=-—(m-l) 
4 À Vi (+ mão O (mão | (4.26) 
Finalmente, passando para o domínio da frequência, vem 
(4.27) 


Wi (f) = At Ri (0) + 2 s 
q= 


= 
Wi (gát) coszrfq4t+ Vi (mAt) cos2xfm st] 
l 
4.2.3. Análogamente ao que foi feito no caso contínuo, vamos determinar a média de con- 


junto destas estimativas. 
(1) Embora estas expressões sejam definidas para |r| <,n, apenas interessa usar valores de r positivos (11% é par) 
desde O até m, em virtude de (4.25). 
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Se fizermos 


v(zãá)=4t 3 o(7—-qát) (4.28) 
q=-— oc 
podemos escrever 


fi () = vy(7,4t) Do (5) D (5) Wi (r4t) 


(4.29) 
= vy (1,4t) D (5) Fr (rá?) 
em que se fez Do Cims)=— eD(tmát)=0 
Tomando as médias de conjunto em (4.29) obtemos 
<U()>=v(gA)D() Y(rAt) 
(4.30) 
= v(,,4t)D (7) F (7) 
Passando para o domínio da frequência, resulta 
<WO>=A (t, To) + O(E) + W(£) (4.31) 
em que 
1 E — j2xf: a q 
A rm = Us dt) é dt = 3 (E — —— : 
( a fo ) T ER “o (4.32) 
Podemos escrever (4.31) de duas maneiras diferentes : 
— Fazendo A (47) “+ W(£) = Wa (f) = s W (4) , temos 
At d=—oo At 
| < Wi) >=0(f) + Wa(f) | (4.33) 


o 
by 


made À (E) OA = Gti 5 Q(t- Th nisi 
At ini Gá At 


| <W()>=Q (E) *W(£) | (4.34) 


A expressão (4.33) aparece-nos como a convolução entre a janela Q (f) e um espectro perió- 
dico Wa (f) resultante duma sobreposição de frequências do espectro a analisar W (f). 

A expressão (4.34) pelo contrário aparece-nos como a convolução entre o espectro pretendido 
W (f) e a janela periódica Qa (f) resultante da repetição (com sobreposição) de Q (f). 


Se W(f) for limitado à frequência E , Wa(f) coincidirá com W'(f) dentro do intervalo 


1 1 4+ ' 4 . . . * 
[E | < pes e a média < W'(f) >> é na maior parte deste intervalo muito aproximadamente igual 


a Q(f) * W(f). Apenas nos extremos os lobos de Q (f) poderão captar contribuições das frequên- 


cias de Wa (f) superiores (em módulo) a a 
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Qa (f) por seu turno só contém um ramo dentro do mesmo intervalo, (excepção feita aos 
valores extremos), coincidindo aproximadamente dentro dele com Q (f). 


, ã fadis PS iai 
No caso porém de o espectro W (f) não ser limitado à frequência DAE dá-se uma sobreposição 
de frequências que pode ser muito importante, desvirtuando mais ou menos profundamente a 
estimativa. 


Este facto, denominado aliasing por Blackman & Tukey, será estudado em pormenor, mais 
adiante. 


4.2.4. Estabilidade 
O cálculo da variância das estimativas, supondo o ruído gaussiano fornece o seguinte resultado 
sen27f n' dt 


Ps : | 
f 1Qa (E — fi) + Qu (E+HE)TW (E rowe-s( a) df df 


di (fi) mai 08 k. E 
| fO (E— 6) Wi) df | 


(4.35) 


|] . * . m 
sendo n = maior inteiro de (n — — 
3 


A expressão (4.35) é susceptível de se simplificar, tal como a sua congénere (4.17), obten- 
do-se para o número de graus de liberdade k: 


— Espectro W (f) pouco acidentado à escala de 


n dt 
És 1 
kma(2 5) t=1,2...m—1 
m 3 
(4.36) 
k = [modo t=0,m 
m 3 


— Espectro W (f) constituído por um pico 
k = 2 


No caso de se realizarem p estimativas, de que se toma a média, há que multiplicar por p os 
valores anteriores. 


4.2.5. Tal como no caso contínuo, apenas nos interessa calcular W' (f) nas frequências 
E : 


Assim, teremos 


à r Wi a úga qrr iii (1) 
Wi ( — = At|Vi(oO)+-2 3 Wi(gAt) cos + Vi (m At) cos r 7 (4.37) 
2m dt q= m 
Esta restrição do domínio de É ao conjunto discreto dos valores — - vai permitir-nos um 


2m 


cálculo simples de Wi re) a partir de Wi (rát). 


2m 


Para isso calculemos W' ( peida usando D (7) = Ds (7). 
2mÃt 


1 A . + a 
(1) Notar que Wi (f) é uma função periódica de periodo —— e simétrica em relação às frequências 
At 
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o ão — —H— re Hc da O TT A A 


wi (5557) = At [vi (0) + 2 Ex ya (q 4t) cos w icdiad + Wi (mAt) cos ee | (4.38) 
2m Â t n q=l m 
pois Vi (q 4t) coincide com W (q At) dentro do intervalo m At. 
Como D (7) Ds (+) =D (7), tem-se 


Wi (rAt) = "Fo (rAt) D (5) 


É (r4t) = To (rAt) D (3) 
e finalmente, no domínio de f, 


Wi (£) = W5 (8) * Q (f) (4.39) 
Esta expressão, quando utilizada para os valores f = É o conduz a uma relação muito 
m 
simples, no caso das janelas Ds» (7), obtendo-se 
1 1 1 
wWy=iwi(y++-wi o) 
2 2 2m At 
- Lai /r=—1 E cad r 1 uai /r+1 
a) = 106 (356) Ev (o) + (5) 
“Aamat 4 A2mdt E Bs 2màdt; cá 4 2m dt (4.40) 
E LM é Sã asõ m—l 


1 Ts fm l.s/ 4 
(a) = (E 
É 24 2 “NYamat Lire “N24t 


resultando para outras janelas expressões formalmente semelhantes. 


: = Gba r a 
Torna-se portanto muito mais fácil calcular wi, Ene) e passar depois para Wi ( a) 
2m dt 2m dt 
segundo as fórmulas anteriores, do que proceder como no caso contínuo, em que se modifica primeiro 
a estimativa da autocorrelação e só depois se aplica a transformada de Fourier. 


4.2.6. A resolução é definida como no caso contínuo, tendo em vista a independência das 
estimativas, resultando 
1 
Af=—— (4.41) 
m dt 
(embora se calculem os pontos r/2m At necessários para definir completamente Wi (f) como já 
frisámos). 


4.3 ORGANIZAÇÃO DO CÁLCULO DO ESPECTRO DE POTÊNCIA 


Temos agora todos os elementos para organizar um projecto de cálculo de W (f), dentro das 
hipóteses em que a análise estatística feita é válida, isto é, no caso dos processos ergódicos, de dis- 
tribuição gaussiana. 

Este projecto, ou melhor, este esboço de projecto será elaborado de acordo com as exigências 
de resolução, estabilidade, duração viável dos registos e ainda largura de banda do processo, no 
caso discreto. 

E dizemos esboço de projecto porque, como veremos, há um certo número de correcções a fazer 
em diferentes alturas da estimativa, correcções essas que são necessárias para obviar à natureza por 
vezes acidentada do espectro e ao inconveniente espaçamento dos dados. 


Assim, e antes de passarmos a estas correcções, podemos estabelecer o projecto da seguinte 
forma; 
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4.3.1. Caso contínuo 


Fixemos a resolução df e a dispersão Dp 
De (4.7) tira-se 


1 
Ta =— 
; | 


De (4.16) e do quadro I obtém-se o valor de k 


=] |— 
D.? 


A 
p 


Admitindo que W (f) não tem variação muito rápida (veremos a correcção adequada mais 
adiante), tem-se, de (4.21) 


Se este valor T' não for demasiado grande para as possibilidades de registo podemos tomar 


um único registo, vindo então 


——— | Tm 


No caso contrário arbitrar-se-á p, obtendo-se para duração de cada registo 


, de (4.19) 


Te = E Ama 
p 3 


Notar porém que, por ser Tm <. Tk se deverá ter 


doar 


Além disso, o uso de T” em vez de T só tem sentido quando Tm << Tk. 


4.3.2. Caso discreto 


Supondo o processo de frequência limitada F (o caso contrário será também visto mais tarde), 
e fixando como no caso anterior Àf e Dp, raciocinaremos como segue: 


Como vimos no $ 4.2.3, convirá ser F = + donde 
t 


A seguir, de (4.41) conclui-se ser 
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-O ORÇAMENTO DUMA EMPREITADA 


A CERTEZA DE 
TEMPOS REDUZIDOS DE PARAGEM 


POR AVARIA DAS MÁQUINAS 
PERMITE-LHE CONSIDERAR CUSTOS MAIS BAIXOS 


o equipamento 


CATERPILLAR 


garante-lhe um máximo tempo de serviço activo 


SOCIEDADE TÉCNICA DE EQUIPAMENTOS E TRACTORES, LDA. 


SEDE PROVISÓRIA: AV. DE MADRID, 3-B-LISSOA 1-TELEFS.: 724053 (PPC A), 710620--TELEGS.: STETRA 
«CATERPILLAR Y KTRAXCAVATOR» SÃO MARCAS REGISTADAS DE CATERPILLAR TRACTOR CO. 


CATERPILLAR 


apresenta 


| 
| 


proni 
0 
TT 
poa 
sm 
MUEGIHIA 


| 
] 
l 


M CAIXA DE VELOCIDADES 
DADA SISTEMA "POWER SHIFT! 


POTÊNCIA CAPACIDADE 


caso VERDE E GUINÉ 


SEDE PROVISÓRIA AV. DE MADRID, 3- B-LISBOA 1 


DISTRIBUIDOR EXCLUSIVO EM PORTUGAL CONTINENTAL, ILHAS AJACENTES, 
SOCIEDADE TÉCNICA DE EQUIPAMENTOS E TRACTORES. L.” 


TELEFS. 724053 (PPC/A)-710620-TELEG. «STETRA> 
* «CATERPILLAR» E «CAT» SÃO MARCAS REGISTADAS DE CATERPILLAR TRACTOR, CO. - PEORIA, ILLINOIS - E.U.A. 


Fazendo intervir (4.36), vem 


pelo que 


sendo n -- 1 o número de pontos por registo. 

De posse de todos os elementos dos $$ 4.3.1 e 4.3.2, não temos mais que efectuar a estima- 
tiva ou estimativas, de acordo com o $ 4.1.1 e com 0 $ 4.2.1 e seguintes. 

No caso de se utilizarem p registos, convém efectuar p estimativas da autocorrelação Wi, e 
com o valor médio obtido calcular Wi. 

O resultado é evidentemente o mesmo que se obtém se determinarmos uma estimativa de Wi 
por registo e efectuarmos a média no fim (em virtude do carácter linear da relação que liga ' a W!) 
mas aquele procedimento tem a vantagem de poupar operações. 


5. O MÉTODO DIRECTO 


54 CASO CONTÍNUO 


5.1.1. Como foi indicado (8$ 2.2.1 e 3.3) multipliquemos agora si. (t) por uma janela B(t), 
função par, e que se anula para |t| > T/2; e efectuemos a estimativa [cf. (2.1)] 


Wi (£) = po (5.1) 


a — jarft 
/ B(t)si(t)e dt 


t 
— T/2 
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que se pode escrever ! 


1 


So 
» — j2rft 
Wi (f) = froso dt 


— OG 


5.1.2. A fim de tomarmos a média de conjunto destas estimativas — provenientes do conjunto 
das funções si (t) — é conveniente servirmo-nos da autocorrelação associada [cf. (2.3)] 


W o à ha à TA E EA 
Peg À B(t- 5) (t-)B(t+5) (t+ ) dt (5.2) 


que está relacionada com W' (f) através da transformada de Fourier. 


Teremos 
Es 1 1 T o T % 
Wi — q— i a A i B RA a 
<< (7) > 3] <s(t De(t+5)> ( ) els h a (5.3) 
ou ainda Pas 
< Wi (1) >= De (1) Y (5) (5.4) 
em que 
Ta 
1 T T 
De ao senta B E fer B = d 
(7) = / (1 =) (t+ -) t (5.5) 
Tp 


De (T) é uma janela equivalente no domínio dos atrasos — no que respeita valores médios, note-se 
— e que se anula para |7|>T. Por isso a relação (5.4) é válida em todo o domínio de 7. 
Passando, como sempre, para o domínio da frequência, de (5.4) obtém-se 


|<WB>=Qe) + w(f) (5.6) 
em que Qe (f), janela espectral equivalente (para determinação de médias) é dada por 


Do 
— j 2rf 7 


Qe (f)=[| De(i)e dz (5.7) 
— oa 
É curioso notar que as expressões (5.5) e (5.7) são formalmente análogas a (2.3) e (2.4). 


É de esperar portanto uma relação correspondente a (2.5). 
Será a seguinte 


1 A — j2rft |2 1 ne 
Qe =] [ B (t) é dt tr AD EE (5.8) 
sendo o 
é — j 2x ft 
1O= [Be dt (5.9) 


(1) É indiferente usarmos como limites de integração + co ou L T/2, em virtude do carácter limitado do domí- 
nio de B(t). 
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5.1.3. Estabilidade 


O cálculo da variância das determinações W (f) conduz (1), no caso gaussiano, ao seguinte 
resultado 


var (W'(fy) |) =<Wi(fi)>º + | [1 E) +) | WE) ae | (5.10) 


mi CIC) 


Resulta assim para a dispersão relativa 47 


* 


= fics) (E + ED | W(£) de —|2 


1 — co 
As(ty)=1 CésadE Paga o. A 
3 (fi) + Tê <Wi(fi) > (5.11) 


Vemos que se tem 4z (fi) > 1, para cada estimativa, qualquer que seja a duração T do registo. 

Pode ver-se que, para janelas B (t) do tipo D (7), a expressão (5.11) conduz a um valor de 47 
que é aproximadamente igual à unidade, para W (f) de andamento gradual, e igual a 2, para fi=0. 

A escolha de B (t) depende principalmente dos requisitos de precisão. 

Com efeito é de esperar que B (t) pouco influa na estabilidade, uma vez que os dados si. (t) 
nos merecem — admitimos por hipótese — igual dose de confiança, ao contrário do que sucede na 
técnica da autocorrelação em que os valores de “Ti. (7) correspondentes aos maiores atrasos são 
calculados com menos dados. 

O problema aqui reside nos lobos laterais introduzidos pelas descontinuidades em t = + T/2 
da janela B, (t). 

Quanto mais gradual for o andamento de B (t) na vizinhança de t==+ T/2, menos oscilações 
introduzirá na sua transformada J (f) e consequentemente em Qe (f). 


5.1.4. Resolução 


Aplicam-se aqui as conclusões tiradas no método da autocorrelação:: 


, : g 1 a 
— Bastará efectuar-se estimativas de W'(f) para valores de f espaçadas de —— , visto que 


WJ! (f) tem como transformada de Fourier a função "Wi (x) que se anula para |7|>>T (ver 5.2). 
— À resolução depende porém da configuração das janelas Qe (f), não se podendo portanto 


exprimir Af de modo genérico em função de a À 


5.22. CASO DISCRETO 


9.2.1. Podemos reduzir este caso ao caso contínuo, precisamente através da janela B (t), 
dando a esta a forma duma sucessão equi-espaçada de funções à (2) (fig. 5-2). 
Ackasu (3) propõe uma janela uniforme idêntica à V (7, dt) D, (7) que utilizámos em $ 4.2.3, 


B (t) = y (t, dt) B, (t) (5.12) 


(1) Ref 2. 

(2) Atitude esta que não era possível no método da autocorrelação, em que apenas as estimativas desta são 
modificadas através de janelas, 

(3) Ref. [2]. 
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Uma vez que em W' (f) só intervêm módulos de transformadas de Fourier, resulta que as nossas 
estimativas não serão afectadas por uma translação no tempo de s' (t) B (t). 

Nestas condições, e por se torrar formalmente mais simples, vamos considerar o intervalo T 
com origem em t = dt. 


8(t) 


IDA pu IM pn , MA 


Fig. 5-2 
Então (5.12.) escreve-se agora 
B (t) = 4t | 3 (48) + s D(t— qÃt + = dlt— (a + 150) (5.13) 
q=? 


Utilizando em seguida a expressão (5.1), e atendendo a que a função st (t) é, com a nova 
origem dos tempos, conhecida desde t= At até t== (n + 1) At, vem para W' (f): 
E du a 
Wi (£) = e [e (E) + b' 9 | (5.14) 
com . 


a (D=4t|+5! (SB cosa rf dt+ > 3 si(gát)cos2r Ffqgát+ Esp Un 40) At) cos 25 (n+1) At | 
q=2 


(5.15) 
[rose [Esp sm semantses s si (qdt)sen 25fq à t+ , kn 4-1) Atlsen ze f(n 41) e] 
q=2 


5.2.2. Resolução 


Tal como no caso contínuo, verifica-se 


<W'(f) > = Qea (f)+*W (f (5.16) 
sendo, para a janela B (t) utilizada 
fAt)sen (nrfÃt) |2 
ea (£) = T q 3 a | 5.17 
Dea (E) | nsen « fÃt 


que é uma função periódica (no domínio de f) de período E « (cf. 4.34) (Fig. 5-3) 
St 


* * . me A * 1 EA 
Vemos assim que vai haver sobreposição de frequências separadas de rr através dos lobos 
t 


principais, desde que 2Fát>>1, em que F é a frequência máxima do processo, suposto de banda 
limitada. (Fig. 5-4) 
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É, de resto, o que se passa, também como já vimos, no método da autocorrelação, e constitui 
uma característica inerente aos métodos de amostragem (sampling). 


A resolução é definida como 


Sf=—— = (5.18) 


tal como no método da autocorrelação, uma vez que se verifica que para este espaçamento as duas 
janelas correspondentes são relativamente independentes. 


Fig. 5-4 


As janelas Q» (f) daquele método, cruzam-se para este intervalo Af a meia altura. 
Aqui teremos aproximadamente o mesmo: 


(1-8 )mom(1= 125) 
2 2ndt 


4 
2 


T (para n>>D1) = 40/ da altura máxima 


, " 
a a ms CE 
n2 2n T 


1 , l A 
Apesar de a resolução ser A há que calcular W! (f) a intervalos , como já indicámos, 


pelo facto de a sua transformada ser de duração limitada. 


5.2.3. Estabilidade 


Supondo não haver a referida sobreposição de frequências (aliasing) o espectro W (f) só con- 
terá um lobo principal. Nestas condições, pode ver-se que a expressão (5.11) da dispersão relativa 
toma os valores mencionados no $ 5.1.3 do caso contínuo, tendo-se portanto, nas zonas em que 
W(£) não apresenta picos, k = 2, valor este que se reduz a metade para f =. 

No entanto, e como já vimos antes, a média de p estimativas diferentes (e independentes) terá 
2p graus de liberdade (p, para f= 0). 
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5.3. ORGANIZAÇÃO DO CÁLCULO DO ESPECTRO DE POTÊNCIA 
5.3.1. Caso contínuo 


Fixemos, como sempre, df e Dp. 
De (5.12) tira-se 


+. 


T=T= 
k AF 


duração de um registo. 


Admitindo dum modo geral k =2 por estimativa, será, como vimos k =2p. O valor k 
tira se de D,, extraído de (4.16), pelo que obtemos o número p necessário de registos 


Temos assim tudo quanto precisamos para o cálculo directo: o número de registos e a dura- 
ção de cada um deles. 

Escolhida a janela B (t), não temos mais do que efectuar as estimativas (5.1) para os p registos 
diferentes e proceder à determinação da sua média aritmética. 


5.3.2. Caso discreto 


Além de 4f e Dp há que estabelecer a frequência máxima em jogo no processo. 
E teremos assim, sucessivamente 


1 
T = Te=— e = — | como no caso anterior. 
A f 2 
E ainda, 
| 1 Es 
At = E espaçamento do sampling necessário 
2 
Tx ; Ki 
a sendo n + 1 o número de pontos exigidos por um registo 
4 t 


N == 


N=(n+1)p<=np| o número total de pontos a utilizar. 


Com estes elementos elaboramos agora as estimativas (5.16) com as quais calculamos W' (f), 


por meio de (5.15), para as frequências múltiplas de a” 


Finalmente, com as p estimativas correspondentes aos p registos independentes necessários, 
determinamos a respectiva média aritmética, que constituirá a nossa estimativa do espectro de potência. 


6. COMPARAÇÃO DOS 2 MÉTODOS 


6.1. Vamos agora representar esquemáticamente os dois métodos estudados, em igualdade 
de resolução 4f, banda F (caso discreto) e estabilidade traduzida pelo número k de graus de liberdade- 
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O quadro II mostra-nos os elementos necessários ao cálculo do espectro, nos 2 métodos, por 
ele se vendo que o número de pontos da função s'(t) a registar, é aproximadamente o mesmo, 
assim como o espaçamento, que só depende da banda do processo. 


QUADRO II 
Dados : F, df, k 


Parâmetro Método da autocorrelação Método da análise directa 
; nnpra k 
N.º de registos p | Arbitrário, sendo 4p << 3k ps 
1 k Geral 
Duração de cada registo Tk Tk = f + E) [AF Tk = 1/Af 
t C 
Atraso máximo Toi Tm = AE e Mc A 
2F 
N.º de atrasos m+1 m = — é 
Af 
1 1 
Espaçamento At At= — At = —— 
2 F 2 F 
Caso 
discreto 
N.º de pontos por registo n + 1 n == A Ê aaa) n= 2F 
Bt 13 2p af 


2F /Pp k 2F kk 
N.º total de pontos N/N=pn=— Es + =) N 
é i af (3 2 


A precisão dos 2 métodos não é exactaménte a mesma, em igualdade de resolução Af. 

Isto resulta de as janelas Qe(f) e Q(f) não serem idênticas. Em geral, para os casos apon- 
tados é Qe (f) pior que Q (f), no que respeita à amplitude dos lobos laterais. 

Quer dizer que para a mesma resolução 4f, as estimativas do método directo vão ser 
mais influenciadas pelos valores do espectro relativos a frequências na vizinhança dos lobos laterais, 
do que no método da autocorrelação. 


A comparação dos dois métodos no que se refere ao custo do cálculo — traduzido através do 
número de operações envolvidas — será tratada noutra publicação. 


6.2. Estudámos comparativamente os dois métodos de determinação do espectro de potência 
W'(f), obtendo-se a maior parte dos resultados (incluindo as normas de projecto) estritamente na 
hipótese de s! (t) constituir ruído gaussiano, e de o espectro não apresentar acidentes muito pronun- 
ciados (à escala da largura das janelas Q (f)). 

Igualmente suposemos, no caso dos registos directos, que a frequência de amostragem era 
dupla duma certa frequência F, acima da qual o espectro se anulava (pelo menos práticamente). 

Pode-se perguntar agora o que sucede se alguma ou algumas daquelas hipóteses se não 
verificarem. 
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No que respeita à banda F ou ao acidentado do espectro, há um certo número de correcções 
que, dentro da aproximação que estes métodos implicam, se podem considerar satisfatórias. Daremos 
adiante uma breve descrição destas correcções. 

Quanto ao facto de o ruído a analisar não ser gaussiano, isso vai traduzir-se por uma maior 
dispersão relativamente aos valores previstos. 

Lembremo-nos, porém, que outras fontes de dispersão existem, como as falhas de estaciona- 
ridade, e que portanto os resultados que atrás indicámos no estudo da estabilidade se devem enten- 
der como aproximados, sendo esta aproximação geralmente satisfatória, segundo mostra a experiência. 


7. OBSERVAÇÕES AO CÁLCULO E CORRECÇÕES 
7.1 CORRECÇÃO PARA OS ESPECTROS DE ANDAMENTO ACIDENTADO. PRÉ-ÊNFASE 


7 1.1. A existência do operador Q (f) « nas nossas estimativas, em que Q (f) 9 (f), implica 
como vimos, contribuição de frequências vizinhas, contribuição essa que pode ser bastante apreciá- 
vel no caso de o espectro W (f) tomar valores muito elevados sobre os lobos laterais das referidas 
janelas (Fig. 7-1). 

Por outro lado, a análise da estabilidade que foi feita mostra-nos que os valores de k usados 
no projecto são válidos apenas na condição de W (f) variar relativamente pouco para um intervalo 
da ordem da distância entre lobos. Caso contrário aquele número k baixa apreciavelmente. 


Por estes dois motivos, uma vez fixada a janela, um caminho a seguir será procurar tornar o 
espectro que se analisa pouco acidentado (1), o que se consegue por meio de filtros, ou, analitica- 
mente, alterando convenientemente os dados. 

Esta alteração dos dados pode realizar-se no caso discreto de modo simples, por meio de 
combinações lineares dos tipos 


s 
Ya= Xq— 5 & Xq-r — combinação linear móvel (7.1) 
p=) 


s 
Ya=X + 5 Er Yq-r —série auto-regressiva (7.2) 
= 


em que Xq representa o elemento de ordem q da sucessão s.i (p At), iniciada arbitrariamente, e Yq O 
elemento correspondente da nova sucessão. 


(1) Embora não se necessite uniformizar absolutamente o espectro, esta técnica (Ref. 1) recebeu a denomina- 
ção de «prewhitening» inspirada na uniformidade do ruído branco (white noise), 
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Vista do inferior de uma das Fábricas "“EFACEC" 
( Arroteia — S. Mamede de Infesta ) 


Distinguem-se nesta fotografia : 
— Um transformador de 150C0 kVA/60 000 V (marca EFACEC) 
com regulação em carga, em curso de fabrico. 


— Um disjuntor de 220000 V (marca SOPREL), aguardando 
ensaio de rigidez. 


— Nota-se também um transformador de ensaios de 400 000 V 


urbar 


imgemieros 


san sebastián 
av generalísimo, 6 
espanha 
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vibradores externos e internos 


réguas vibrantes equipamentos vibrantes 


para superfícies 


energia electrica, peneumatica e gasolina 


alta frequência e meia frequência 


representante 


sociedade comercial romar ltda 
rua da boavista, 83, 1º d 


lisboa 


“+ 
as 
q 
> 
-Q 
3 
a 


Prova-se facilmente que se tem (1) 


q s 12% 2 
W'a (E) =| q 5 a jsfrat (7.3) 
Wa (£) r=1 
Wa (f) s —jzmfrat|—2 
=|1—35 E 7.4 
ig zé | (7.4) 


respectivamente, em que Wa (f) e W'a (f) são os espectros periódicos (v. $ 4.2.3) correspondentes 
às sucessões Xq e Ya. 
Estes dois tipos de combinações lineares, quando devidamente associados conduzem a expres- 


a a cl di , 1 
sões racionais em cos27=fAt de qualquer grau, e portanto periódicas em f com período a 
t 


Claro que uma vez feita a estimativa baseada nos novos valores há que multiplicar o espectro 
obtido pelo factor inverso do utilizado. 


7.1.2 Quando o espectro apresente linhas, isto é, picos de amplitude muito grande comparada 
1 
com o restante espectro, e de largura muito pequena comparada com a banda o problema de obter 


uma combinação linear conveniente pode tornar-se inviável, em virtude do grande número do pará- 
metros 2, e É, necessários. 

Nestas condições torna-se preferível efectuar duas análises dos dados, uma estimativa na zona 
que contém o pico, e outra longe desta. 

Esta segunda determinação implica uma pré-ênfase sem grande rigor, que nos permite analisar 
Wf) no restante domínio, sem que o pico vá causar perturbação. 

Quanto à avaliação de W(f) no próprio pico, uma simples estimativa comparada com a que 
referimos pode dar uma boa ideia da amplitude e largura do pico. 

Para melhor localização do pico pode-se, evidentemente, usar um filtro de banda pelo qual se 
faz passar o espectro a analisar. 


7.1.3. Quando os picos do espectro aparecem nas muito baixas frequências, há um tipo de 
correcção muito mais simples. É este o caso das variações lentas (no tempo) e termos constantes 
introduzidos durante as medidas e que os sistemas digitais não cortam. 

A fim de obviar a este inconveniente, e porque não conhecemos o valor médio e a baixa 
frequência (estranhos ao processo) introduzidos, Blackman & Tukey (2) aconselham, no método da 
autocorrelação, a subtrair sistemáticamente, dos valores | (q 47) calculados, expressões convenientes. 

Assim, por exemplo, a rejeição do valor médio pode ser obtida subtraindo de cada estimativa 
Vi (g At), o quadrado da média aritmética dos n + 1 valores Xq = s.i(q At), ou seja 


1 n+1 —a 
Rr >, si(a dt) | (7.5) 


Esta subtracção, vindo afectar igualmente todos os valores de 'Vi(q At), irá afectar a estimativa 


Wi (f) apenas nas frequências f=0 e f= a (para as janelas D (7) usadas). 


(1) Estas relações só são exactas se Xq for definido para q == — co. No caso real, se Xg se iniciar p. ex. para 
q=-0,a relação 7.1 só tem sentido para Yq a partir de q==s. Para q >s a não utilização de alguns elementos vai in- 
troduzir em 7.3 soluções não periódicas. 

(2) Ref. [1] — págs. 47. 
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Prova-se que, para o ruído branco, isto equivale estatisticamente a multiplicar o espectro do 
processo que se pretende estudar, isto é, o processo sem o valor médio introduzido, pelos factores 


(1-5) (1-5) (7.6) 


respectivamente, nas frequências f=-0 e f=1/2 m Àt pelo que há que multiplicar as determinações 


assim feitas pelos factores inversos daqueles. 
Pode-se, alternativamente, subtrair de cada observação si (qdt) a média de todas as observações. 


ç] n41 
—— E 84 (gát) (7.7) 

com as mesmas correcções a fazer. 
Para rejeitar simultâneamente valores médios são nulos (termo constante) e variações lentas, 


indica-se a subtracção às estimativas ''r (q At) da expressão 


PE -1 2 


Es 3 E 3 
: NERDS | e q 
Sp + ST. fls] (7.8) 
16 « / 
em que si. ; - e si. representam as médias aritméticas respectivamente dos n + 1 valores 


si (q dt), do terço superior e do terço inferior da mesma sucessão. 
Para esta correcção os factores finais a usar serão os inversos de 


“ 
x 


(1- E) pisa) nas seg pf = (7.9) 
n n 2m dt 


7.2. EFEITOS DUM ESPAÇAMENTO INCORRECTO St 
7.2.1. Vimos, quando estudámos o método da autocorrelação para o caso discreto, que se tinha 


|<WO>=0 (8 * Wa (b (4.33) 


Ao estudarmos o método directo obtivemos, no caso discreto, 


<W'(f >= Qe (8) *W (£) 


em que Qen (f) era uma função periódica de período Y simétrica em torno de f = = : 
Escrevendo 
Qu () = A(£ 1) * Qo(f = 30 ( fa E 
obtemos uma expressão idêntica a (4.33) 
| <Wif)> = Qe(f) * Wa (f) | (7.10) 


Portanto, em qualquer dos métodos, a média das estimativas pode apresentar-se como um 
integral de convolução entre uma janela Q (f) ou Qs. (f) não periódica e um espectro periódico, 
Wa (f) resultante de W (f) por sobreposição de frequências. 
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Para melhor compreendermos o significado destas expressões, suponhamos o registo s.i (t) 
ilimitado. Nestas condições, tanto D (7) como B(t) se podem considerar de duração ilimitada. 


Usando janelas uniformes, teríamos em qualquer dos casos Q (f) = Qe(f) = 9(f). 
As expressões (4.33) e (7.10) escrever-se-iam então 
< Wi(f) > = Wa (f) (7.11) 


Quer isto dizer que, mesmo nesta hipótese duma amostragem indefinida (n =) há inde- 
terminação para W (f). 

O que se determinaria univocamente (em média) seria apenas a versão Wa (f) do espectro 
pretendido W (f). 

Wa (f) representa portanto, no limite, tudo quanto podemos aspirar a determinar. 

Esta indeterminação resulta da natureza discreta dos nossos registos. 

Sendo as funções s! (t) apenas conhecidas (ou utilizadas) a intervalos iguais At, sabemos da 


teoria da interpolação que este conhecimento é suficiente para determinar todos os pontos intermé- 


dios (e daí W (£f)) apenas no caso de o espectro da função (e com ele, W (f)) ser limitado a F = k.. 


2Ãt 
1 
Se for este o caso, Wa (f) coincide com W (f) dentro do intervalo |f| < CA e a expres- 
são (7.11) escreve-se 


<Wi()>=W(f) (7.12) 


que é a relação ideal (cf. caso continuo com T= o), 


E arca E duras 1 
No caso porém em que W (f) não é limitado, ou é limitado, mas F + CÊ teremos que 
t 


— ou nós conhecemos F> ——, e o espectro obtido (no limite) Wa (f) é uma sobreposição 
t 


dos valores do espectro pretendido W (f) em diferentes frequências 
— ou nós não conhecemos F, e W (f) então é completamente indeterminado, não podendo 
nós sequer fazer uma ideia da validade da nossa estimativa. 


7.2.2. Uma interpretação geométrica curiosa deste espectro Wa (f) (1) resulta de se escrever 
esta função sob a forma 
Wa (t)=WI(t)+HW (2fn— t)+W(QEn+t)+W(s fx-f)+H+W(4 fx +f)+.. (7:13) 


visto W (f) ser uma função par, sendo 


1 
fx = (7.14) 
2At 


Se dispusermos os vários argumentos no domínio de f que se irão sobrepor em Wa (f), obte- 
remos para uma frequência f fx (o caso f>>fwx é redutível a este) a disposição da Fig. 7-2. 


ca 
6 
WN sy 
4f N 
4 
Yo he Yo f 21 

O O O N ” 

H f 

I 

Fig. 7-2 Fig. 7-3 


Se dobrarmos em zig-zag este eixo das frequências pelos pontos q É x, obteremos agora a dis- 
posição da Fig. 7-3, em que se vêem sobrepostas as frequências com que W (£) irá contribuir para Wa (f). 
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Wa (£) aparece-nos assim como resultante da dobragem do espectro pelas frequências múltiplas 
de fx (folding frequency ou frequência de Nyquist) (1) (Fig. 7-4). 


Fig. 7-4 


7.2.3. Em resumo: Analiticamente ou geométricamente vemos como, no caso de o espectro 
W (f) ser limitado a F = fx, todos os termos do desenvolvimento de Wa (fi) são nulos, à excepção 
do primeiro, que é precisamente W (fi), sendo portanto o folding inoperante. Nestas condições a 
média de todas as estimativas possíveis, que é igual a Q(f) * Wa (f), só recebe contribuição das 
frequências dentro da resolução Af, e apenas para fi próximo do extremo fw, é afectada pela 
periodicidade de Wa (£) (Fig. 7-5). 


E ia! 7 % 


Fig. 7-5 


No caso de F> fx, porém, a média referida pode ser completamente estranha ao pretendido. 
Bastará para tanto que algum dos termos W (2q fx + fi) tenha valor apreciável (v. Fig. 7-4). 


7.3. PROBLEMAS DE ECONOMIA DO CÁLCULO. AS BAIXAS FREQUÊNCIAS 


7.3.1. A indeterminação de W (f) que acabamos de estudar constitui, por ser inerente aos 
processos de amostragem, um inconveniente que, para uma dada banda F e uma vez escolhido 
o espaçamento Àt, é impossível evitar. 

Vamos mostrar, no entanto, que em determinadas condições se pode evitar o aumento daquela 
indeterminação, decorrente duma redução de frequência de dobragem fx. 

Isso consegue-se, como veremos, quando se está interessado apenas na zona inferior do espectro 
e se pretende melhorar aí a resolução, sem onerar demasiado o cálculo (2), 

Antes porém vamos estabelecer um método de, a partir dum dado registo discreto, fazer variar 
os factores estabilidade, resolução e número de multiplicações necessárias, à custa duma diminuição 
da frequência fx. 


(1) Ref. 11). 

(2) O volume do cálculo, quando este é feito por meio de computadores digitais, depende largamente do número 
de multiplicações necessárias, podendo dizer-se que este número, que é dado aproximadamente por M=mn por 
registo, no método de autocorrelação (Ref. 1, pág. 44) é uma boa medida do custo do cálculo. 
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O referido método consiste no seguinte: 
Tomemos um registo de n +41 dados espaçados de à4t e desdobremo-lo em h registos de 
n+1 


dados espaçados de 47 ==h At, como se indica na figura. 


Fig. 7-6 


Apliquemos o método da autocorrelação e comparemos uma estimativa do tipo (a) que utiliza 
todos os pontos e um atraso máximo m dt: 


1) com 1 estimativa a partir de um dos h registos do tipo (b) que utiliza apenas pontos 


de hát em hát, com um atraso máximo m' À 7 
2) com a média das h estimativas possíveis idênticas à anterior. 


Em vista de análise feita, pode obter-se o quadro seguinte. 


QUADRO III 


1 estimativa (a) 


Estabilidade | Resolução Custo Dobragem 
ko | MF M fx 
1 estimativa (b) h estimativas (b) 
Estabilidade Resolução Custo | Estabilidade | Resolução | Custo 
k' = k Af=Af |IM=Mh'|m=m/h|k'=kh |Af'==Af | M' =M/h 


K=kh |Af=Afh |M=M/h m=m k=k [Af=4fh IM =M 


e | À À E | TT 


k=k/ht AP=Agh? M=M im =mh | K = k/h |Af=4f/h?|M' =Mh 


Para todos os casos: fx ==fn/h 


Finalmente, determinemos a relação entre o espectro periódico Wa (f) correspondente ao registo 
primitivo, e o espectro Wan (f) correspondente a qualquer dos registos (b). 
Tem-se 


Wii = d (6 o) * W(£) 
At 


Wan (£) = An (6 E + W(£f) 
Mt 


A (6 E) =) (€- E = 33 (E—29£n) 


q=-oo 


em que 
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/ à Ma 2 1 
An(£ ac) = 2 A (E as fx) 
hãt, t 


+ 


E portanto 
Vialli É Vs ( + dE fx) (7.15) 
r=o0 


O espectro periódico Wan (f) resulta assim de h sobreposições do espectro original Wa (f). 


har Ap (8, Vhat) 
f 
f 
+ valsas) 
f 
Fig. 7-7 


7.3.2. Vamos agora descrever a compensação do abaixamento de fx, nas baixas frequências, 


referida no último parágrafo. 
Suponhamos que pretendemos levar o cálculo de W (f) até uma frequência mais baixa do que 


2mdt 
tendemos baixar este, mantendo embora a resolução. Desde que estejamos interessados só nas 
baixas frequências podemos utilizar o método do parágrafo anterior. 

Consultando o quadro III vemos que, mantendo a estabilidade podemos obter o que preten- 
demos com uma estimativa e m'== m/h, ou com h estimativas e m'==m. No 1.º caso mantemos a 
resolução e baixa-se o número de operações de M para M/hº e no 2.º baixa-se a resolução, man- 
tendo-se M. 

Em qualquer dos casos porém f'x = fx/h. 

Se o valor fx utilizado já for o menor que não prejudica as determinações no terço inferior 
do espectro, ou se F- 3fx/2, então não podemos pensar em reduzir fx (por ex. a metade) sem 
tomar qualquer precaução. 

Essa precaução consiste no seguinte: 

Efectuemos a combinação linear 


Za = Xa - Xq-t + .... + Xq-k ti (7.16) 


sendo Xq o elemento da ordem q da sucessão s.. (qi). 
Tem-se como se sabe 
Wa (8) RA a fantast 


1 da 
Wa (f) o 


1 , 
aquela que o nosso projecto de cálculo determinou (a) , sem aumentar o cálculo, ou que pre- 


2 


)= Y (£) (7.17) 


sen rfkdt 


sen nfdt 
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em que Wa (f) diz respeito à sucessão s.. do processo em estudo, e W'a à sucessão definida por 7.16. 
Apliquemos agora a combinação (7.16) ao espectro representado na Fig. 7-8. 
Vemos que Y (f) atenua considerávelmente o espectro Wa (f) na região em que ele não coin- 
<ide com W (f) (dentro do intervalo [o, + f=], claro) (Fig. 7-9). 


Fig. 7-9 Fig. 7-10 


Se agora utilizarmos apenas 1 em cada 2 dados Z,, obteremos um novo espectro periódico 
Wa (E) = Wa (8) + W'a (E + fx) (v. 7.15) cuja frequência de Nyquist fx é igual a metade da 
anterior f'N = fx (Fig. 7-10). 

Se finalmente multiplicarmos W“a (f) 1) por 1/Y (f), vemos que o resultado obtido difere 
muito pouco do pretendido Wa (f) dentro da zona |f| < fx/4, crescendo a divergência entre os 
Wa (£) 

Y(£) 
Portanto dentro da zona |f| < fw/1 (ou um pouco mais além) os efeitos da diminuição de fx pelo 
uso de metade dos valores Z, foram muito atenuados, o que não nos deve surpreender, pois lem- 
bremo-nos que estes Zy mesmo reduzidos a metade contêm quase toda a informação das baixas frequên- 
cias, O que não aconteceria se nos limitássemos a utilizar metade dos Xq (Fig. 7-10), caso em que se 
obtém W, muito diferente de Wa. 
Em resumo, querendo economizar o cálculo devemos proceder como segue: 


valores e Wu(f) rápidamente, à medida que |[f| varia de fn/4 até fx. (1) 


1) Efectuar a combinação (7.16). 
2) Aproveitar estes dados, utilizando um sim, um não. 


(1) Wa 


ia Wa . sec? 7 fAt, como se verifica facilmente. 
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3) Realizar uma estimativa de W' seguindo os métodos estudados. 
4) Multiplicar a estimativa final pelo inverso de Y (f). (1) 


Com uma estimativa obtemos assim a mesma estabilidade e resolução, e, o que é mais importante, 
reduzimos o número de multiplicações de 4 vezes (as operações suplementares provenientes do 
uso de Y(f) pouco afectam esta economia). 

O contra está na perda de informação nas frequências mais elevadas. 

Com as h estimativas 'neste caso, 2) e um número de atrasos duplo (h vezes maior), conse- 
guimos descer na resolução para metade (sempre nas baixas frequências), mantendo a estabilidade 
e não onerando o cálculo. 

O método é generalizável (2) a k>8 (k = 4h), sendo o 1.º lobo obtido, ou seja, para 


Ele = , O menos afectado. 


8. ESTIMATIVAS PRÉVIAS 


8 1. Vistos dum modo geral os tratamentss que os registos devem sofrer, antes de se efec- 
tuarem as estimativas, e as respectivas correcções finais, resta dizer que toda esta análise pressupõe 
em geral um conhecimento aproximado prévio do espectro a estudar, a fim de que se possa proceder 
à pré-ênfase e estabelecer a resolução requerida. 

Esse conhecimento pode obter-se por uma medida ou por uma estimativa grosseira, existindo 
métodos de cálculo numérico aproximado baseados em simples somas e diferenças dos valores 
«discretos dados. 


8.2. Um exemplo destas estimativas grosseiras (estimativa piloto) é o indicado no quadro IV (3), 

Como vemos, parte-se dos dados todos, cujo espaçamento se supõe adequado à banda F. 
-e calcula-se a soma dos quadrados das diferenças de valores contíguos. 

Esta quantidade depende largamente do intervalo (fx , fx/2). 

Passa-se para a oitava seguinte (fn/2, fx/4) pela mesma técnica, somando préviamente os 
“valores 2 a 2 e tomando apenas 1 sim e 1 não, o que equivale a duplicar o espaçamento. 

E assim sucessivamente até ficarmos reduzidos ao quadrado da soma de todos os dados — estima- 
tiva da banda mais baixa. 

Os valores obtidos são (muito grosseiramente) proporcionais às potências relativas aos sucessivos 
“intervalos de oitava, a partir de fx no sentido descendente (fig. 8-1). 

Para melhores aproximações existem métodos mais precisos (4), 


9, CONCLUSÕES 


As normas gerais que acabamos de descrever não podem constituir um critério exacto de avaliação 
-de todo e qualquer espectro de potência. 

Não constituem um critério exacto, porque como já referimos, existem causas de erro que 
-derivam do facto real de os processos em estudo se afastarem mais ou menos dos modelos matemá- 
“ticos que lhes atribuimos. 
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(1) Note-se que se tem < Wi(f)> =Q f) « Wa(f)e que multiplicar Wi (f) por 1/Y(f) não é o mesmo que 
“fazê-lo para Wa (f) a menos que Y(f; seja constante. Como porém o domínio da convolução é praticamente muito 
pequeno comparado com fN/4 , Y (rt) pode tomar-se como constante em cada banda, o que valida o método. 

(2) Ref. [1], pág. 44. 

(3) Ref. [1], pág. 46. 

(4) Ref. [1], pág. 135. 
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Este facto deve ser bem presente sempre que procuremos interpretar os resultados obtidos. 

Por outro lado não se pode estabelecer um critério universal a menos que nos confinemos aos 
espectros pouco acidentados (à escala da frequência maior a analisar) e de frequência limitada bem 
definida F, para os quais a análise feita parece ser bastante satisfatória. 

Na hipótese contrária, ou até por considerações de outra ordem que só cada caso pode impor, há 
que estabelecer, por exemplo, uma frequência máxima a analisar, desdobrar a análise em várias estima- 
tivas de diferente resolução e até estabilidade diversa, conjugar o estudo analítico com o experimental, 
enfim, um sem número de soluções que a prática e a finalidade de cada caso possam sugerir. 
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A II CONFERÊNCIA MUNDIAL DE ENGENHARIA 
SÍSMICA 


por J. FERRY BORGES 
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Chefe do Serviço de Edifícios e Pontes do 
Laboratório Nacional de Engenharia Civil 


O presente trabalho foi apresentado durante a «Reunião sobre Construção Anti-Sís- 
mica» efectuada no Laboratório Nacional de Engenharia Civil em 27 e 28 de Fevereiro 
de 1961. Nele se relatam as orientações seguidas e os principais resultados obtidos pela 
II Conferência Mundial de Engenharia Sísmica, no estudo dos seguintes problemas : 

a) Comportamento das estruturas durante os sismos e equipamento para o estudo de tal 


comportamento. 


b) Movimentos vibratórios do solo durante um sismo. Influência dos tipos de terrenos 


de fundação. 


c) Dimensionamento de construções anti-sísmicas. Regulamentos. Sismos recentes e danos 


consequentes, 


1 — INTRODUÇÃO 


A II Conferência de Engenharia Sísmica teve 
lugar em Tóquio em Julho de 1960 (1) e deu con- 
tinuidade a uma primeira conferência realizada 
na Califórnia em 1956. 

Participaram nesta segunda conferência cerca 
de 600 engenheiros e outros especialistas, tendo 
sido apresentados perto de 120 trabalhos. Estes 
trabalhos tiveram qualidade e originalidade tal 
que bem pode afirmar-se que a conferência mar- 
cou o início duma nova época neste ramo de en- 
genharia. 

A aplicação de conceitos muito recentes e de 
novos instrumentos de investigação permitiu 
obter conclusões que terão reflexos técnicos de 
grande importância. De facto, os fenómenos sís- 
micos sômente são convenientemente explicados 
recorrendo à mecânica das vibrações aleatórias 
que tem sido desenvolvida nos últimos anos. 
Sômente, também, os computadores analógicos 
e digitais constituem a ferramenta adequada à 
análise teórica do comportamento dinâmico de 
estruturas complexas. 

Finalmente, únicamente graças a aparelhos de 
realização também muito recente é possível efec- 
tuar estudos experimentais satisfatórios. Com- 
preende-se pois que num futuro próximo se as- 
sista a um grande desenvolvimento da engenha- 
ria sísmica. 


Graças a uma comparticipação da Fundação 
Calouste Gulbenkian, que cumpre agradecer, teve 
o Autor possibilidade de participar na Conferên- 
cia, para a qual preparou uma comunicação (2) 
em que se discutiam as orientações mais conve- 
nientes para o estabelecimento da regulamenta- 
ção de construção anti-sísmica. Durante a Con- 
ferência apresentou também resultados quantita- 
tivos relativos aos coeficientes sísmicos de vá- 
rias construções danificadas durante o sismo de 
Agadir, resultados esses de interesse para o jul- 
gamento dos valores dos coeficientes sísmicos 
correntemente adoptados (3). 

O grande trabalho dispendido pelo Japão, para 
a organização da conferência, traduziu-se em 
enorme eficiência e proveito para os participan- 
tes. Além das sessões de trabalho, que se desen- 
volveram simultâneamente em duas salas, teve 
lugar uma exposição, principalmente de equipa- 
mento, e efectuaram-se ainda visitas aos vários 
institutos de investigação de Tóquio em que se 
estudam problemas sísmicos. 

Reuniões efectuadas entre delegados dos dife- 
rentes países conduziram à demonstração da ne- 
cessidade de constituir uma associação interna- 
cional que permitisse dar continuidade a estas 
reuniões. Tal associação está sendo formada e, 
graças a ela, será possível efectuar a próxima 
conferência mundial. 

O estudo que se apresenta, além de procurar 
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sintetizar os principais resultados obtidos na 
conferência referida, pretende fornecer um pano- 
rama geral da evolução que têm sofrido os prin- 
cipais problemas de engenharia sísmica. 


2 — COMPORTAMENTO DAS CONSTRU- 
ÇÕES DURANTE UM SISMO E EQUI- 
PAMENTO PARA O ESTUDO DE TAL 
COMPORTAMENTO. 


A matéria correspondente a este assunto pode 
ainda subdividir-se conforme se trata de análise 
teórica (em regime elástico ou plástico) ou de 
análise experimental (incluindo estudos sobre 
modelos e protótipos) e instrumentos para o 
estudo do comportamento dinâmico. Analisa-se 
portanto seguidamente, cada um destes problemas. 


2.1 — Análise teórica do comportamento das estrutu- 
ras em regimen elástico 


Os esquemas teóricos simples assimilam as 
estruturas a osciladores elásticos, com ou sem 
amortecimento, e de um ou vários graus de 
liberdade. 

São numerosos os estudos recentes que con- 
sideram o comportamento dum oscilador harmó- 


estudos iniciais de Housner (4) que bem se pode 
considerar terem praticamente constituído a base 
da engenharia sísmica moderna. 

Para descrever o comportamento dum oscila- 
dor, interessa introduzir a noção de espectro. 

Considere-se que se dispõe do registo do mo- 
vimento do terreno, numa dada direcção, obtido 
por um sismógrafo, fig. 1a). A partir de registos 
de acelerações pode passar-se para velocidades 
ou deslocamentos por conveniente integração. 
Analise-se a resposta dos osciladores harmóni- 
cos amortecidos de diferentes períodos próprios 
e com diferentes graus de amortecimento, quando 
o terreno sofre acelerações de acordo com o dia- 
grama indicado. Designando por z (t) a acele- 
ração do terreno, é fácil demonstrar que a velo- 
cidade relativa do oscilador em relação ao ter- 
reno pode ser calculada pelo integral também 
indicado na fig. 1. Tal integral permite calcular 
os valores máximos dos deslocamentos, veloci- 
dades ou acelerações relativas. Por conveniência, 
considerem-se as velocidades. Marcadas, num 
diagrama, estas velocidades máximas em função 
do período, obtém-se um conjunto de pontos, 
que uniremos por segmentos de recta. O dia- 
grama assim obtido designa-se por espectro da 
velocidade máxima, fig. 2. 


ONA MS 
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cos wit—<) dz 


Sv = max X = espectro de velocidade 
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Z 
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3 — fracção do amortecimento crítico 


Sd = — Svy — espectro de deslocamento 


2" 
2x 


Sa = FT Sv — espectro de aceleração 


Fig. 1 — Comportamento dum oscilador harmónico 


nico, de um grau de liberdade, com diferentes 
graus de amortecimento, sob a acção de movi- 
mentos vibratórios registados quando de sismos 
reais, devendo mencionar-se, em particular, os 
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Este espectro indica, pois, quais as velocida- 
des máximas sofridas pelas estruturas simples 
do tipo considerado. Interessa frisar que, na gama 
de períodos compreendidos entre 0,5 e 3 segun- 


